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Doublon-Holon Bound States in a Half-filled 
Band 

In a recent Letter, Sakai, et al.Q] performed a numer- 
ical study of the Hubbard model on a square lattice in 
attempt to command a clear view of a doped Mott insula- 
tor. Among the results they present as being new is the 
" non-trivial finding" that the spectral weight immedi- 
ately above the chemical potential in the lower Hubbard 
band, A(x), increases rapidly with hole doping (x) at the 
expense of spectral weight in the upper Hubbard band. 
Although this result is presented without reference to any 
prior work, it is far from new. It was first observed in a 
cuprate by Chen, et al.[2jin 1991 and was immediately 
explained by Eskes, et al.[3j as a ubiquitous phenomenon 
in a doped Mott insulator. In fact, the density of states 
in Fig. (3a) of Sakai, et al.[l[ is identical to that in Fig. 
2 in the second of Ref. 3. 

Sakai, et al.pj also propose that A(x) is described by 
2x + 2nd (counting both spins) where rid is the number 
of doubly occupied sites. While the dynamical part of 
A(x) is certainly mediated by double occupancy, it is 
not equivalent to it as can be seen from the perturbative 
expression, 


A(x) = 2x 
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derived by Harris and Langejj] in 1967 (also not cited). 
This expression contains a static part (2x) because each 
hole creates two addition states at low energy and a dy- 
namical part which is strictly positive and proportional 
to the expectation value of the kinetic energy. In fact, 
all orders of perturbation theory^, 0] increase A(x) be- 
yond its atomic limit of 2x. The result of Sakai, et al.[lj 
follows from the perturbative analysis only if the virial 
theorem holds for the Hubbard model, which it does not. 
As a consequence, any relationship between the dynam- 
ical part of A(x) and double occupancy is not one of 
equivalence but rather supervenience. 

In fact, Sakai, et al.[l[ are well aware that their anal- 
ysis is problematic. They note correctly, that on their 
account, A(x = 0) = 2nd at half-filling, which is non-zero 
rather than vanishing as would be required by the Mott 
gap. To correct this flaw in their interpretation of A(x), 
they state, without any kind of proof, that some type of 
bound state between doublons and holes saves the day 
by transferring spectral weight back to the upper Hub- 
bard band. It is unclear what this statement means since 
there is no energy scale associated with nd- It is simply 
a number. Further, Sakai, et al.[l[ offer no explanation 
as to why their mechanism for the pseudogap, screening- 
induced degradation of the binding energy, would lead 
to an abrupt rather than the more physically reasonable 


gradual reduction of the Mott gap. 

However, the claim that bound states between dou- 
blons and holons mediate the Mott gap is not without 
merit because without them, the holons (arising entirely 
from quantum fluctuations) would be mobile thereby de- 
stroying the Mott gap and preventing antiferromagnetic 
order. As a consequence, establishing that such bound 
states are the propagating degrees of freedom in a half- 
filled band is the essence of the Mott problem, provided, 
of course, that the band structure of these elementary ex- 
citations yields the gapped spectrum indicative of Mot- 
tness. To this end, the Wilsonian program of integrat- 
ing out the high energy scale indispensable and has been 
recently carried out exactly Q for a half-filled band de- 
scribed by the Hubbard model. The exact low-energy 
Lagrangian demonstrates that while there are no bare 
degrees of freedom which propagate, composite excita- 
tions corresponding to bound doublon/holon or double 
hole/electron pairs do propagate. The spectral weig htg 
of the doublon/holon pairs turns on for U/2 — At < u> < 
U/2 + At while that for the double- hole/electron pair is 
non-zero for —U/2 - At < cj < -U/2 + At, thereby de- 
scribing the dispersions for the upper and lower Hubbard 
bands, respectively. For U < 8t (8t the bandwidth) the 
gap for the composite excitations (not that for the elec- 
trons) vanishes giving rise to a discontinuous jump in the 
density of states at the chemical potential to the corre- 
sponding value in the non-interacting system. At finite 
doping, such bound states still persist and generate a 
pseudogap r 5j , though with a new energy scale not related 
to the Mott gapvia the screening mechanism envisioned 
by Sakai, et al.[lj. 
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